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CHAPTER  1 
INTRODUCTION 

Periodic  tuning  is  a  virtually  essential  aspect  in  the  maintenance 
of  almost  any  physical  plant.  One  method  of  such  tuning  is  that  known  as 
sensitivity  points  turning  method.  The  theory  of  this  method  and  application 
to  single-input/single-output  (SISO)  systems  have  been  fairly  well  developed. 
This  study investigates  the  application  of  sensitivity  points  turning  concepts 
to  the  tuning  of  multi-input/multi-output  (MIMO)  systems. 

The  concept  of  sensitivities  in  feedback  systems  was  originally 
introduced  by  Bode  in  the  1940 's  [1]^  Methodologies  for  using  sensitivity 
information  for  optimal  tuning  of  system  parameters  date  back  to  the  late 

1950 's  and  early  1960 's  and  include  work  done  by  Meissinger  [2,3]^  Margolis 

/  : , 

and  Leondes  [4,5],  and  numerous  others  (see  [6,7]).  These  works  all  involve 
the  use  of  system  sensitivities  in  order  to  provide  a  searchless  means  of 
computing  the  gradient  of  a  given  cost  function  which  possesses  some 
extremal  value  that  is  related  in  some  way  to  some  desired  characteristic 
performance  of  the  system  under  tuning. 

Included  in  this  class  of  optimization  techniques  are  the  iterative 
parameter  optimization  techniques,  such  as  those  of  Meissinger  [2,3],  Brunner 
[8],  Roberts  [9],  and  Kokotovic  [10].  The  choice  and  generation  of  pertinent 
sensitivities  for  use  in  such  techniques  are  topics  treated  in  works  by, 
among  others,  Kokotovic  [10],  and  Wilkie  and  Perkins  [11,12].  The  sensitivity 
points  tuning  process,  developed  by  Kokotovic  [10],  is  an  iterative  process 
which  obtains  sensitivities  in  a  straightforward  manner  and  then  adjusts 
parameters  in  an  optimizing  fashion. 


1.1.  Sensitivity  Points  Tuning  Concepts 

Consider  the  system  depicted  in  Figure  1.  In  this  arrangement, 
which  occurs  frequently  in  practical  applications  of  control  systems, 
compensator  blocks  K^(cu,s)  and  ^(B^.s)  of  known  structure  are  used  for 
control  of  system  blocks  W^(s)  and  W^Cs) ,  whose  structures  need  not  be 
necessarily  known. 

For  a  given  input  r(s)  and  some  small  deviations  Aa^  and  AB^ 
from  some  nominal  parameter  values  a?  and  8?,  respectively,  a  Taylor  series 
expansion  of  the  output  response  yields 

y(a°  +  Aou  ,  0°  +  ABj  ,  s) 


=  y(aj  ,  8°  ,  s)  +  Z  ^ 

1=1  l 


9/0 

—  y(a  ,  s) Aa  + 


+  38.  ’  s^A0i  +  - • • 

3=1  J  J  J 


where 


9  A 

- — y(a.,s)  =  sensitivity  of  the  output  y(a.,B.,s)  with  respect 
°  i  1  3 


the  parameter  ou  at  =  a^,  and 


-T-—  y(8.,s)  =  sensitivity  of  the  output  y(ot.,8.,s)  with  respect 
3D  J  1  J 

J  0 
to  the  parameter  8.  at  3.  =  S.. 

J  J  J 

Note  that  (1)  may  be  rewritten  as 
+  Aa.  ,  8°  +  AS .  ,  s) 


,r(a.  +  Act.  ,  3.  4-  AS.  9  s) 
i  l  j  J 

■  y<»°  .  .  ■)  *  V<“?  •  •  s)  '  + 


y (a?  ,  3°  ,  s)  •  A3  +  ... 


If  the  expansion  is  truncated  at  the  first-order  term  and 
y(ou,8j,s)  is  taken  to  be  the  nominal  response  of  the  system,  then  the 
derivation  of  the  system  response  yfc^  +  Aa^.B^  +  AB^  ,s)  from  the  nominal 
response  y(a?,8^,s),  denoted  as  the  error  e(s),  may  be  approximated  by 

e (s)  38  y(a?  ,  6°  ,  s)  •  Aa 
1  1 


+  v^y(o‘°  ,  8*?  ,  s)  •  AS 
c  1  J 


(3) 


If  a'  =  a^  4  Aa,  3'  =  8^  +  AS,  and  Aa  is  known  to  be  small,  then 


the  approximation  in  (3)  may  be  rewritten  as 


e(s)  38  V  y (a T  ,  ST  ,  s)  •  Aa 
a  i  j 


4  V  y(a'  ,  S'  ,  s)  •  A8 
B  1  1 


(4) 


The  approximation  made  in  (4)  implies  that  knowledge  of  the 
sensitivities 


r— —  y(aT,s)  and  y(3',s)  would  allow  one  to  approximate 
da.  1  a  B  J  J 


appropriate  values  of  Aa^  and  AS^ ,  subtract  them  from  the  values  of  a^ 
and  3'  resident  in  the  system,  and  reduce  the  amount  of  error  e(s).  The 
resultant  e(s)  would  have  a  magnitude  approximately  that  of  the  higher- 


order  terms  ignored  in  the  approximation  of  e(s)  made  in  (4).  The  magnitude 


5 


of  this  resultant  e(s)  would  be,  therefore,  smaller  than  that  of  the  original 
e(s).  Repetition  of  this  process  through  several  iterations  should  reduce 
the  magnitude  of  e(s)  to  some  minimal  value  (the  presence  of  measurement 
noise  may  make  it  impossible  to  reduce  the  magnitude  of  e(s)  to  zero) . 

Note  that,  due  to  the  Laplace  integral's  property  of  differentiability 
with  respect  to  system  parameters,  this  process  is  also  directly  applicable 
in  the  time  domain. 

The  sensitivities  required  for  this  procedure  may  be  obtained  in 
the  following  manner.  Refer  again  to  the  system  diagrammed  in  Figure  1. 

The  system  transfer  function  and  sensitivities  are 


y(or  ,  3  , s)  = 


W1(s)K1(ai,s) 

1  +  W1(s)K1(a. ,s)W2(s)K2(S. ,s) 


W  (s)K  (a:,s)W  (s)K  (S's) 

~  v(a'  s)  - - - - — - - - - j  •  vT  K  (a;,s)  •  r(s) 

i  1  (1  +  W1(  J)K1(a;T,s)W2(s)K2(3j  ,s))4  3“i 

W1(S)^T  Ki(Vs) 

1  /  X 


1  +  W1(s)K1(a^,s)W2(s)K2(6jr5s) 


.  3a.  K1(VS) 

_ 1 _ #  l _ 

1  +  W1(s)K1(a',s)W2(s)K2(£;,s)  K^otT.s) 


y  (  jl'  ,  3 ' ,  s ) 


1  +  W1(s)K1(a',s)W2(s)K2(2j,s)  ’  la.  ln  Kl(ai,S> 


•y(a^,3j,s) 


98  y(6j»s) 

j 


W1(s)K1(a',s)W2(s) 


(1  +  W1(s)K1(a',s)W2(s)K2(8j,s))  j 


2*  3TK2(6rs)  °r(s) 


W  (s)K  (ar,s)  3 

1  ,w1(3)Kl(a-,s)w2(s)K2(5]r,irw2(s>'?ir  K2<ej.s)y«.;,6-,s) 


(8) 


W1(s)K1(a',s) 


•  W0(s)K„(B;,s) 


1  +  W1(s)K1(a',s)W2(s)K2(8',s)  2  ^2 v 

3 


3«7K2<“i’,) 


W1(s)K1(a£,s) 


T  •  W.(s)K,(B;,s) 


1  +  W1(s)K1(a^,s)W2(s)K2(Bjr,s)  2v°'~2VMj 


jt *"  V6rs)  • 


(9) 


The  reason  for  using  the  sensitivity  filters  — —  £n  K. (a.,s)  and 

e  7  3a .  1  l 

i 


3  3 

T-r—  £n  K„(37,s)  in  (7)  and  (9),  respectively,  instead  of  -r —  K,  (a',s)  and 

OD  .  Z  1  oOL  .  1  1 

]  i 


t-t—  K  (3r,s),  as  in  (6)  and  (8),  is  to  avoid  using  higher-order  sensitivity 

0  O  ,  /  1 

1 


filters  when  the  parameter  is  in  the  denominator  of  the  compensator  transfer 
function . 

Example:  Let 


f  (s) 

V8i’s)  = 


(10) 


where  f  (s)  and  f  (3C,s)  are  polynomials  in  the  Laplace  operator  describing 
2  P  ] 

the  zeros  and  poles,  respectively,  of  K^fS^s),  Then 


where 


ith  (6rs)  = 

j 


f  (S) 
z 


f  (s:,s) 


(fP(6j 's)>2  p  3 


(ii) 


tt-  S-n  K  (e.",s) 
1  J 


1  _1 .  -  V  (O'  s) 

K„(B;,s)  38,  j*  ' 


2  j 

f_  (s) 


J 


ytpr  TT 


(12) 


The  order  of  the  filter  in  (11)  is  lower  than  that  in  (10)  by  the  order  of 


The  signals  described  in  (7)  and  (9)  may  be  obtained  as  diagrammed 
in  Figure  2.  The  points  and  are  known  as  the  "sensitivity  points"  of 
the  system  at  which  appropriate  signals  may  be  easily  picked  off  for  use  in 
calculating  sensitivites . 

1.2.  Parameter  Deviation  Estimation 

The  approximations  of  the  parameter  deviations  Lou  and  Aii  may  be 
obtained  in  any  number  of  ways  in  order  to  satisfy  whatever  error  minimization 
criteria  that  are  specified.  One  of  the  more  common  tuning  goals  is  that  of 
finding  and  A£L  to  minimize  the  magnitude  of  difference  between  the 

system  response  and  a  nominal  reference  response  over  the  entire  response 
period  T0 .  For  this  purpose,  one  may  use  a  verv  common  cost  function;  namelv. 


SENSITIVITY  FILTER 


Configuration  for  obtaining  sensitivities  of  the  system 
diagrammed  in  Figure  1. 


that  of  rms  error  over  the  specified  response  period.  The  cost  function 


may  be  expressed  as 


J  =  /  (e(t)Tdt 
0 


(13. 


In  discrete  time,  (11)  becomes 


N 


J  =  Z  (e(kT))i 
k=0 


(14 


where  NT  =  T  . 

R 


Now  let 


a 

"Act- 

Y  = 

and  Ay  = 

_8_ 

„AS_ 

If  a  time-domain  analogy  is  made  of  the  error  approximation 
expressed  in  (4) ,  the  tuning  goal  in  each  iteration  becomes  that  of  finding 


N 


Ay*  =  arg  min  [  E  (e(kT)  -  V  y(a . , 6 . ,kT) *Ay)  ] 
Ay  k=0  Y  i  J 


(15) 


in  discrete  time,  or 


Ay*  =  arg  min  [  /  (e(t)  -  V  v(a . , S . t) • Ay )  dt ) 

Ay  0  r  1  2 


(16) 


in  continuous  time. 

If  data  processing  were  to  be  performed  by  a  digital  computer  or 
microprocessor,  the  tuning  objective  in  each  iteration  would  be  that  stated 
in  (15),  which  may  be  interpreted  as  a  minimization  problem  in  a  Euclidean 
vector  space  of  order  N+l.  In  this  context,  (15)  may  be  written  as 


*  •  J 


J 


where 


Ay  =  arg  min  J SAy  - y 
Ay 


vYy(a.,8.,0) 


S  =  V  y (a . , 0 . kT) 
Y  1  J 


_AYy(ait8j ,NT)_ 


U  =  e(kT) 


e(NT) 


The  solution  for  (17)  may  be  given  by 

Ay*  =  (sWy  (18) 

T  —1 

if  (S  S)  exists  (implying  that  no  linear  dependencies  exist  between  the 
columns  of  S)  [13].  Note  that  for  n  +  m<2,  (18)  always  applies;  for 
n  +  m>  3,  the  possibility  of  linear  dependencies  between  columns  of  S  arises 
and  a  more  complex  computation  than  (18)  must  be  undertaken  in  order  to 
ensure  that  a  suitable  solution  to  (17)  may  be  obtained. 


1.3.  Notes 


Up  to  this  point,  this  discussion  has  mentioned  only  the  use  of  a 
first-order  method  of  parameter  approximation.  Higher-order  methods  for 
approximation  may  be  used  [14],  but  generally  do  not  yield  significantly 


superior  performance.  In  addition,  higher-order  methods  usually  require 
longer  computation  times  in  each  iteration. 

A  second,  more  important  point  to  be  made  here  is  that  the 
sensitivity  model  used  to  obtain  the  pertinent  sensitivities  need  not  be 
the  system  itself  or  a  copy  of  the  system,  as  might  be  inferred  by  (7)  and 
(9),  and  as  depicted  in  Figure  2.  If  a  full  model  of  the  system  is 
available  (in  either  digital  or  analog  form),  this  model  may  be  used  as 
the  sensitivity  model.  In  addition,  since  the  parameter  deviations  of 
interest  are  only  being  approximated  at  each  iteration  (as  opposed  to 
one  single  exact  calculation  for  the  entire  process) ,  a  sufficiently 
descriptive  reduced-order  model  may  be  substituted  as  a  sensitivity  model 
if  the  full-order  model  of  the  system  is  inappropriately  large  or  too 
complex  for  such  use.  If  the  reduced-order  model  is  well-chosen,  the 
performance  of  the  sensitivity  points  tuning  process  will  not  be  significantly 
degraded  by  such  approximation.  If  the  chosen  reduced-order  model  is  suspect, 
however,  two  changes  may  be  made  to  the  procedure  outlined  above  in  order 
to  reduce  effects  of  using  such  a  mode. 

1.  The  sensitivity  point  may  be  located  in  the  original  system,  as 
depicted  in  Figure  3,  instead  of  the  reduced-order  model,  whereupon  only  the 
input  and  output  of  the  reduced-order  sensitivity  model  need  be  accessible. 
This  eliminates  the  need  for  the  internal  structure  of  the  sensitivity  model 
to  be  identical  to  that  of  the  system.  Matching  only  of  input-output 
character istics  of  the  sensitivity  model  and  system  would,  therefore,  be 
sufficient  in  formulating  a  viable  sensitivity  model. 


Alternate  configuration  for  obtaining  sensitivities  of 
the  system  diagrammed  in  Figure  1. 


*V-|*  I 


2.  The  estimation  of  the  parameter  deviation  vector  Ay  may  be 
made  using  a  steepest  descent  estimate  of  the  form 


*  T 
Ay  =  cS  p 


c  s  m 


instead  of  the  least-squares  solution  given  in  (18).  This  update  form  is 
less  sensitive  to  inaccuracies  in  the  sensitivity  matrix  S.  In  addition,  if 
c  is  appropriately  chosen,  this  steepest  descent  update  will  tend  to  have 
faster  initial  convergence.  One  disadvantage  of  this  form  is  that,  with 
constant  c,  the  updates  may  oscillate  about  the  optimal  solution.  One 
possible  solution  to  this  problem  is  to  utilize  the  steepest  descent 
estimate  for  several  initial  iterations,  and  then  use  the  least-squares 
estimate  when  in  the  neighborhood  of  the  optimal  solution  in  order  to  ensure 
convergence. 


1.4.  Subjects  for  Discussion 

In  the  ensuing  discussion,  the  concepts  of  the  sensitivity  points 

a 

method  will  be  extended  to  apply  to  multiple-input/multiple-output  (M1M0) 
systems.  The  Mimo  sensitivity  points  method  will  then  be  used  to  examine, 
first  of  all,  techniques  of  tuning  MIMO  systems  in  a  rather  general  sense, 
and,  then,  the  applicability  of  the  sensitivity  points  method  to  decentralized 
tuning  of  MIMO  systems.  Examples  and  simulations  will  be  included  to  help 
illustrate  points  to  be  made.  Finally,  an  example  of  a  realized  system  will 
be  given  and  discussed. 


CHAPTER  2 


SENSITIVITY  POINTS  TUNING  IN  MULTI-INPUT/MULTI-OUTPUT  SYSTEMS 

In  this  chapter,  the  sensitivity  points  tuning  concepts  described 
in  Chapter  1  for  single-input/single-output  (SISO)  systems  will  be  extended 
to  apply  to  multiple-input/multiple-output  (MIMO)  systems.  Techniques  for 
sensitivity  points  tuning  of  MIMO  systems  will  then  be  discussed. 


2.1.  MIMO  Sensitivity  Points 

The  deviations  for  the  sensitivities  of  MIMO  systems  progress  in  a 
manner  quite  similar  to  that  described  in  Chapter  1  for  SISO  systems.  Consider 
the  very  general  MIMO  system  depicted  in  block  diagram  form  in  Figure  4. 


This  system  has  an  input  vector  R(s)  of  dimension  dR;  output  vector  y(a,8,s) 


of  dimension  d^;  system  feedforward  block  W^(oi,s)  with  transfer  function 


matrix  of  dimension  d  x  d  ,  and  with  adjustable  parameter  vector 

y  K 


a  =  [a, , . . . ,a  ]  ;  and  system  feedback  block  W_(3,s)  with  transfer  function 
In  Z 


matrix  of  dimension  dR  x  d^  and  with  adjustable  parameter  vector 


6  =  3  ]  .  The  transfer  function  of  the  system  is 

I  m 


,-l. 


y(o,3,s)  =  [I  +  WL(a,s)W2(3,a)]  WL(a,s)R(s) 


(19) 


If  one  were  to  follow  a  notation  convention  analogous  to  that  used 
in  the  SISO  sensitivity  derivation,  one  would  find  the  corresponding  MIMO 


sensitivities  to  be 


jf-  y(a',6-,s)  -  -  [I  +  W  (<*-,s)W  (B'.s)]  ^(ct'.s)  •  g|-  W2(S',s) 
j  j 

•  [I  +  W1(a',s)W2(8%s)]“1W1(a%s)R(s) 

=  -  [I  +  W1(a>,s)W2(8',s)]'1W1(a%s)  W2(B\s)y(a,B%s)  . 

j 

(21) 

A  realization  for  obtaining  the  signals  described  by  (20)  and 
(21)  is  diagrammed  in  Figure  5.  If  one  were  to  refer  back  to  the  SISO 
realization  depicted  in  Figure  2,  one  would  see  that  a  commutation  of  the 
sensitivity  filter  and  sensitivity  model  in  Figure  2  would  result  in  a 
degenerate  form  of  the  realization  shown  in  Figure  5 . 

The  lack  of  commutativity  in  the  MIMO  case  has  one  major  implication; 


namely,  that  only  one  sensitivity  filter  ?t  a  time  may  be  used  in  the 
sensitivity-generating  process.  In  the  SISO  case,  commutativity  of  terms 
allows  for  the  placement  of  multiple  filters  on  the  output  of  the 


sensitivity  model.  Sensitivities  with  respect  to  several  different  parameters 


may,  therefore,  be  obtained  simultaneously.  Since  no  such  commutativity 
exists  in  the  MIMO  case,  however,  and  since  only  one  input  vector  signal 
may  be  applied  to  the  sensitivity  model  input  at  any  time,  only  the 
sensitivities  with  respect  to  one  parameter  will  be  generated  at  any  given 
time.  This  parameter  is  the  parameter  corresponding  to  the  sensitivity 
filter  in  use.  If  the  sensitivities  with  respect  to  p  parameters  are  desired, 
then  the  process  of  inputting  y(s)  (or  E(s) ,  depending  upon  which  is 
applicable)  into  a  sensitivity  filter-sensitivity  model  pair  must  be 
executed  p  times. 

If  one  were  to  use  the  MIMO  system  depicted  in  Figure  6,  one 
would  find  the  system  transfer  function  to  be 

y(a,3,s)  -  [I  +  W1(s)K1(a,s)W2(s)K2(6,s)]“1W1(s)K1(a,s)R(s)  .  (22) 


The  corresponding  system  sensitivities  would  be 


-f-  y(a\3',s)  =  -  [I  +  W1(s)K1(a',s)W2(s)K2(S',s)]  1 
i 

•  W1(s)  K1(a",s)W2(s)K2(B',s) 

i 

•  [I  +  W1(s)K1(a',s)W2(s)K2(6%s)]“1W1(s)K1(ct%s)R(s) 

+  [I  +  W1(s)K1(a',s)W2(s)K2(S',s)]'1W1(s)  K  (a',s)R(s) 

i 

=  [I  +  W,  (s)K.  (a',s)WJs)L(S',s)]'V  (s)  -f-  K,(a',s)E(s) 


E(s)  =  R(s)  -  W2  (s)K2  (B'"ss)y  (a  ' ,  8 "  ,s) 

and 

rj~  y(o,8,s)  =  -  [I  +  W1(s)K1(a,s)W2(s)K7(6,s)]“1 

j 

•  W1(s)K1(a,s)W2(s)  ~~  K2(S,s) 

•  [I  +  W1(s)K1(a,s)W2(s)K2(B,s)]_1W1(s)K1(C1,s)R(s) 

=  -  [I  +  W1(s)K1(a,s)W2(s)K2(B,s)]“1W1(s)K1(a,s) 

W2(S)  IT  K2(6’s)y(a,8’s)  •  (24) 

j 

A  realization  for  obtaining  the  signals  described  by  (23)  and 
(24)  is  diagrammed  in  Figure  7 .  An  alternate  realization  of  (23)  and  (24) 
is  shown  in  Figure  8.  The  attractiveness  of  this  alternate  realization 
lies  in  the  fact  that  the  system  itself  may  be  used  as  the  sensitivity  nodel 
without  any  modifications  such  as  the  addition  of  auxiliary  inputs  shown  in 
the  realization  of  Figure  7. 

In  addition,  if  a  reduced-order  model  were  to  be  used  as  the 
sensitivity  model,  it  would  only  need  to  match  the  input/output 
characteristics  of  the  system. 

Use  of  the  alternate  realization,  however,  is  limited  to  those 
svstems  where  K^(a,s)  exists  and  is  realizable,  and  where  W^(s)  is 
readily  replicable.  In  either  case,  as  in  the  more  general  realization 
shown  in  Figure  5,  a  lack  of  commutativity  implies  that  onlv  the  sensitiviti 


K(«\3) 


Configuration  for  obtaining  sensitivities  of  the  system 
diagrammed  in  Figure  6. 


te  configuration  for  obtaining  sensitivities 
diagrammed  in  Figure  6 . 
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with  respect  to  one  parameter  may  be  obtained  at  any  one  time,  since  only 
a  single  sensitivity  filter  may  be  used  in  series  with  the  sensitivity 
model.  One  possible  solution  to  this  problem  would  be  to  use  multiple 
sensitivity  models  (or  multiple  approximate  sensitivity  models),  which 
would  allow  for  the  extraction  of  sensitivities  with  respect  to  as  many 
parameters  as  the  number  of  sensitivity  models  used. 


2.2.  The  MIMO  Sensitivity  Filter 

At  this  point,  it  may  be  appropriate  to  discuss  the  interpretation 


of  the  MIMO  sensitivity  filter, 


9(0 


K(*,s).  In  the  S1S0  case,  realization 


of  the  sensitivity  filter  amounted  to,  quite  simply,  the  realization  of  an 
SISO  transfer  function.  In  the  MIMO  case,  however,  the  transfer  function  of 
the  sensitivity  filter  is  a  matrix  K(-,s)  whose  elements  mav  be  denoted  as 
K_(',s)  (i  denoting  row  assignment,  j  denoting  column  assignment); 

— yy  K(*,s)  is  simply  a  matrix  of  the  same  dimension  as  K(*,s)  whose 


i-j-th  element  is 
Example :  Let 


!(•)  ij 


K. . (• ,s) 


K  =  i 


then 


1 

s+: 


k(-:,s)  = 


o 

o 


o  : 

i 

°! 

0  l 


ith  row  of  K(  • ,  s)  and  use  this  sununed  signal  as  the  ith  element  of  the 

o(  •) 


With  S  and  p  defined  in  this  fashion,  Ay  may  be  estimated  with 


the  same  methods  used  in  the  SISO  case. 

It  should  be  noted  at  this  point  that  none  of  the  outputs  need 
necessarily  be  tuned  with  respect  to  all  parameters.  If,  for  instance,  one 
wished  to  tune  output  y^  by  using  all  parameters  in  y  except  for  some  6  , 
all  that  one  would  need  to  do  would  be  to  set  ~~  y.(t)  =  0  for  all 

OP.  1 

3  3 

t=  [0,Tn]  .  This  would  set  -r-r—  y.(t)Ag.  =  0,  meaning  that  it  would  have 

K  d  p .  1  1 

j 

no  effect  in  the  tuning  process.  The  estimation  would,  therefore,  ignore 

* 


in  approximating  Ay  . 


2.4.  MIMO  Tuning  Techniques 

If  one  were  to  assume  that  the  MIMO  system  to  be  tuned  was,  in 
fact,  a  network  of  SISO  or  MIMO  subsystems,  then  tuning  techniques  would 
fall  into  two  general  categories:  centralized  tuning  and  decentralized 
tuning . 

2.4.1.  Centralized  tuning 

In  this  scenario,  subsystems  of  the  overall  system  are  interconnected 
by  inter-subsystem  feedback  as  well  as  any  coupling  inherent  in  the  plant.  As 
depicted  in  Figure  9,  the  main  purpose  of  a  centralized  tuner  would  be  to 
maintain  the  inter-subsystem  feedback  elements  in  a  way  that  would  assure 
overall  system  stability.  These  particular  feedback  elements  have  less  effect 
on  subsystem  dynamics  than  local  feedback  elements,  but  may  become  important 


if  inherent  coupling  is  strong.  Otherwise,  a  perturbing  input  to  one 


SUBSYSTEM 


subsystem  may  excite  resonances  in  some  less-damped  subsystems  via  inherent 
coupling,  which  may  then  result  in  loss  of  overall  system  stability. 

Example :  Magnetic  supension  system.  Consider  the  magnetic  suspension 
system  sketched  in  Figure  10.  The  objective  is  to  maintain  the  bar  at  a 
nominal  distance  ho  and  to  keep  it  level.  The  linearized  state  equations 
of  the  system  are 


o 


and  b  =  7-^ 

rt 

o 

Notice  that  the  system  may  be  treated  as  two  weakly  coupled  sub¬ 
systems.  Note  also  that  these  subsystems  appear  in  the  system  matrix  as 
diagonal  blocks,  which  allows  for  a  feedback  controller  matrix  of  the  form 


In  this  controller  form,  those  elements  in  the  diagonal  blocks  are  associated 
with  local  feedback  and  the  off-diagonal  elements  govern  cross-feedback. 

Assume  now  that  a=l,  m=2,  g=l,  h  =1,  L=2,  and  that  the  sampling  interval 
is  50  msec.  Assume  further,  for  the  sake  of  simplicity,  that  k^^  =  k2^  =  0.9, 
^12  =  k24  =  ’  and  ^4  =  ^22  =  ^ ’  so  we  will  be  tuning  with  just  the  two  cross 

feedback  elements  k^  and  This  describes  an  unstable  fourth-order 

system.  Computer  simulation  will  demonstrate  that,  for  a  step  disturbance 
applied  to  subsystem  1  with  ^^3  =  ^21  =  ’  t^s  system  is  unstable.  Tuning 

k^  and  k?1  by  the  sensitivity  points  method  with  a  second-order  reference 
step  response  (poles  at  s  =  -2±j3.46)  results  in  convergence  of  k^  and  k^^ 
to  stability  in  12  iterations.  Responses  change  from  unstable  to  32% 
overshoot  in  y^  and  150%  overshoot  in  y^  (tuning  inputs:  unit  step  into 
subsystem  1  (left),  half  step  into  subsystem  2  (right)).  The  tuning  history 
of  this  particular  run  is  given  in  the  Appendix  as  Table  I.  The  effectiveness 
of  a  centralized  tuner  is,  therefore,  quite  clear  and  may  be  critical,  as 
demonstrated  by  this  simple  example. 


2.4.2.  Decentralized  tunin 


Consider  now  the  problem  of  decentralized  tuning  at  each  subsystem 
In  contrast  to  centralized  tuning,  local  tuners  may  use  only  locally 
available  signals  for  the  tuning  process,  as  depicted  in  Figure  11.  This 
implies  that  the  task  of  local  tuning  is  limited  to  the  tuning  of  only  local 


compensators  (which  may  be  in  either  the  feedforward  or  feedback  path) . 
Consequently,  only  those  parameters  in  the  local  compensators  need  be 
tuned . 

One  major  issue  in  decentralized  tuning  is  whether  or  not  exact 
tuning  is  possible  from  the  local  level.  If  the  full  model  of  the  overall 
system  is  available  to  the  subsystem  as  a  sensitivity  model,  such  tuning 
is  possible. 

In  order  to  visualize  why  this  is  true,  assume  that  the  transfer 
function  matrix  of  the  system  may  be  written  in  a  form  such  that  each 
subsystem  lies  within  a  block  on  the  main  diagonal.  This  would  result  in 
the  placement  of  local  compensator  terms  in  blocks  of  the  main  diagonal 
of  their  respective  transfer  function  matrices,  which  would  correspond  to 
the  location  of  their  respective  subsystem  blocks  in  the  main  system 
transfer  function  matrix.  The  structure  of  the  MIMO  sensitivity  filter 
(discussed  in  Section  2.2)  would,  then,  ensure  that  only  local  signals 
would  be  needed  for  local  tuning,  irrespective  of  the  level  of  coupling  to 
other  subsystems  and  irrespective  of  inputs  to  other  subsystems.  This 
implies  that,  if  the  full  model  of  the  system  is  available  as  the  sensitivity 
model,  no  approximation  of  sensitivities  will  be  made;  the  exact  sensitivities 
described  in  (20)  and  (21)  (or  (23)  and  (24))  will  be  generated  instead. 
Availability  of  such  exact  information  ensures  that  exact  tuning  from  the 
local  level  is,  indeed,  possible. 

Irrespective  of  whether  or  not  a  full  system  model  is  available  as 
a  sensitivity  model,  decentralized  tuning  may  be  performed  under  any  of 
several  schemes.  All  of  such  schemes,  however,  seem  to  be  one  or  the  other 
or  a  combination  of  two  basic  methods. 


Tables  II  and  III  of  the  Appendix. 


CHAPTER  3 


DECENTRALIZED  TUNING  OF  MIMO  SYSTEMS 
WITH  SIMPLIFIED  SENSITIVITY  MODELS 

Given  an  MIMO  system,  consisting  of  coupled  subsystems,  consider 
again  the  problem  of  decentralized  tuning  of  individual  subsystem  controller 
sections.  As  depicted  in  Figure  11,  local  tuners  have  only  locally  available 
signals  for  tuning  information. 

It  has  been  established  in  Chapter  2  that,  using  the  MIMO  sensitivity 
points  tuning  procedure,  exact  local  tuning  at  a  particular  subsystem  may  be 
obtained  irrespective  of  coupling  or  of  the  size  and  timing  of  inputs  to 
other  subsystems  if_  a  full  model  of  the  overall  plant  is  available  at  the 
subsystem  for  use  as  a  sensitivity  model.  Such  a  full  model,  however,  may 
not  be  readily  available  for  such  use.  The  problem  to  be  addressed,  then, 
is  whether  or  not  a  subsystem  can  be  satisfactorily  tuned  if  only  a  greatly 
simplified  model  of  the  overall  system  is  available  for  use  as  a  sensitivity 
model . 

In  order  to  assess  limitations  of  decentralized  tuning  with 

a 

simplified  sensitivity  models,  three  tuning  situations  have  been  examined. 

1.  Tuning  with  the  rest  of  the  plant  tuned.  In  this  case,  the 
subsystem  of  interest  is  tuned  up  from  some  detuned  state,  while  all  of  the 
other  subsystems  in  the  plant  are  operating  at  their  respective  fully  tuned 
states.  Simulation  will  show  that,  with  this  particular  situation,  tuning 
with  a  simplified  sensitivity  model  is  as  effective  as  tuning  with  a  full 


model  of  the  plant  as  the  sensitivity  model. 
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2.  Alternate  tuning.  This  situation  involves  cooperation  between 
subsystems.  Subsystems  are  tuned  one  after  the  other  in  a  tuning  "run" 
through  the  entire  plant.  This  process  is  then  repeated  until  successive 
tuning  runs  yield  no  change  in  the  respective  parameters  being  tuned  in  each 
subsystem  controller.  The  only  limitation  of  this  tuning  process  is  that, 
while  a  particular  subsystem  is  being  tuned  in  the  midst  of  a  tuning  run,  the 
rest  of  the  plant  must  be  stable. 

3.  Tuning  with  the  rest  of  the  plant  mistuned.  This  situation,  like 
Situation  1,  does  not  involve  cooperation  between  subsystems.  For  the  purposes 
of  this  thesis,  the  term  "mistuned"  will  have  the  connotation  that  some 
subsystem  not  being  tuned  has  been  detuned  to  such  an  extent  that  the  subsystem 
being  tuned  tunes  to  parameter  values  different  from  those  arrived  at  when 
using  a  full  model  of  the  plant  for  the  sensitivity  model.  This  situation 
arises  when  one  or  both  of  two  conditions  are  met: 

(a)  The  rest  of  the  plant  is  unstable;  or, 

(b)  The  rest  of  the  plant  has  a  natural  undamped  frequency  below 
some  low  frequency  which  is  specific  to  each  overall  plant  and  may  be 
specific  ot  each  subsystem. 

Example :  Magnetic  suspension  system.  Consider,  again  the  magnetic  suspension 

described  in  Chapter  2.  The  controller  matrix  is,  again. 
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As  before,  those  elements  in  the  diagonal  blocks  are  associated  with  local 
feedback  and  the  off-diagonal  elements  govern  cross-feedback.  Since  we  are 
dealing  with  decentralized  tuning,  assume  that  the  off-diagonal  terms  are 
zeros.  For  the  values  of  the  suspension  model  parameters,  assume  that  a=l, 
m=2,  g  =  l,  hQ  =  l,  L  =  2 ,  and  that  the  sampling  interval  is  50  msec.  Tuning 
kn»  ^23’  anc*  k24  the  sens^t^vi-ty  points  method  with  the  full  plant 

model  as  the  sensitivity  model  and  with  a  second-order  reference  step 
response  (double  poles  at  s=-3.0)  results  in  convergence  of  parameters  to 
kn  =  k23=1.16  and  k^^  =  =  0 . 68  .  Assume  for  each  subsystem  that  its 

respective  simplified  sensitivity  model  is  defined  by  those  portions  of  the 
plant's  system,  control,  and  output  matrices  which  pertain  to  that  particular 
subsystem  and  local  feedback  about  that  subsystem.  For  the  left  subsystem 
(Subsystem  1),  then,  the  sensitivity  model  is  defined  by 


12 J  |_2b  '  ^kl!  “kllJ  |_712 J 


S  =  [1  0]  fj. 


The  right  subsystem  (Subsystem  2)  sensitivity  model  is  similarly  defined. 

Detuning  Subsystem  1  alone  to  k^  =  0.5  and  k  =  1.5  and  tuning 
result  in  a  demonstration  of  Situation  1  (tuning  with  the  rest  of  the  plant 


tuned) .  The  sensitivity  points  tuning  procedure  yields  convergence  in  4 


iterations  to  k^=1.15  and  ^2  =  0-69.  The  tuning  history  for  this  tune-up 
is  shown  in  the  Appendix  as  Table  IV. 

Detuning  Subsystem  1  to  k^  =  0.5  and  k^2  =  l*5,  detuning  Subsystem  2 
to  k22  =  0»25  and  k2^  =  0.1,  and  tuning  up  the  two  subsystems  in  alternate 
fashion  result  in  a  demonstration  of  the  alternate  tuning  scenario  of 
Situation  2.  The  sensitivity  points  tuning  procedure  yields  covergence  in 
16  iterations  in  2  "runs"  to  k^=1.2,  k^2  =  0-70,  ^23*^'^*  k2^  =  0.66.  The 
tuning  history  for  this  tune-up  is  shown  in  the  Appendix  as  Table  V. 

Consider  now  the  case  of  tuning  with  the  rest  of  the  plant  mistuned 
(Situation  3).  For  this  example,  let  us  tune  Subsystem  1  while  Subsystem  2 
is  tuned  to  some  arbitrary  state  in  order  to  observe  what  Subsystem  2  must 
be  tuned  to  for  it  to  be  considered  mistuned.  The  poles  of  "the  rest  of  the 
plant"  in  this  case  are  those  of  Subsystem  2,  and  may  be  calculated  from  the 
sensitivity  model  system  matrix  of  Subsystem  2  (in  this  particular  example). 
Now  let  us  assume  the  following  about  all  of  the  tune-ups  associated  with 
this  particular  situation: 

1.  The  tuning  reference  input  to  Subsystem  1  is  a  unit  step. 

2.  The  reference  response  for  Subsystem  1  is  a  unit  step  response 
of  a  second-order  system  with  poles  at  s=-3. 

3.  The  disturbance  input  to  Subsystem  2  is  0.5  step  delayed  by 

0 . 6  sec  . 

A  summary  table  of  initial  and  final  values  of  k^  and  k^9, 
corresponding  values  of  k,.j  and  k7^,  and  corresponding  resultant  undamped 
natural  frequencies  of  Subsystem  1  is  given  in  the  Appendix  as  Table  VI. 
Respective  tuninc  histories  are  shown  in  the  Appendix  as  Table  VII.  As 
can  be  seen  in  Tibi?  VI,  undamped  natural  frequencies  in  Subsystem  2  of  2 


or  above  result  in  k  and  k^  converging  upon  values  in  the  neighborhood 
of  those  attained  when  tuning  with  a  full  model  of  the  plant  as  the 
sensitivity  model,  while  undamped  natural  frequencies  of  1.414  and  below 
result  in  final  values  which  are  much  different. 


CHAPTER  4 


CONCLUSION 

A  general  framework,  has  been  presented  for  the  tuning  of  MIMO 
systems  by  use  of  the  MIMO  sensitivity  points  method.  Examination  and 
computer  simulation  of  the  MIMO  sensitivity  points  method  have  revealed 
that,  for  a  system  composed  of  weakly  coupled  subsystems,  decentralized 
tuning  may  be  satisfactorily  performed  even  if  only  simplified  sensitivity 
models  are  available  for  use.  In  addition,  it  has  been  established  that, 
if  the  full  system  model  is  available  for  use  as  a  sensitivity  model, 
exact  tuning  of  a  subsystem  from  the  local  level  is  possible,  irrespective 
of  coupling  to  other  subsystems  or  inputs  to  other  subsystems.  One  natural 
extension  of  these  results  which  may  be  of  interest  for  further  study  is 
that  of  the  application  of  MIMO  sensitivity  points  principles  for  truly 
adaptive  (instead  of  periodic)  on-line  tuning. 


APPENDIX 


TUNING  HISTORIES  OF  TUNING  EXAMPLES 


TABLE  I.  CENTRALIZED  TUNING  EXAMPLE:  TUNING  HISTORY 


FOR  THIS  RUN:  SUBSYSTEM  1  INPUT  ■  UNIT  STEP; 

SUBSYSTEM  2  INPUT  -  0.3  STEP. 

CROSS-FEEDBACK  GAINSt  K21  AND  K13 

ITERATION 

I  K21  K13 


r  . 


o 

- 1 . 00000 

-1.00000 

l 

-0.8 1344 

-0.79887 

2 

-0.61201 

-0.57941 

-0.3 2783 

-0.23119 

4 

-0. 07316 

+0.29346 

5 

-0.04473 

+2.06399 

6 

-0. 04061 

+2.31908 

7 

-0.03874 

+2.49874 

S 

-0.03907 

+2.35895 

9 

-0.03870 

+2.53242 

10 

-0. 03889 

+2.55015 

1 1 

-0.03877 

+2.53984 

12 

-O. 03884 

*2.54607 

TABLE  II.  DECENTRALIZED  TUNING  EXAMPLE:  ALTERNATING 
TUNE-UP  TUNING  HISTORY 


LOCAL  FEEDBACK  GAINS: 

Kll  AND  K12  CORRESPOND  to  SUBSYSTEM  1 
K23  AND  K24  CORRESPOND  TO  SUBSYSTEM  2 

ITERATION 

1  Kll  K12  K23  K24 

FOP  This  RUN:  REFERENCE  INPUT  -  UNIT  STEP: 

DISTURBANCE  INPUT  ■  0.3  STEP  APPLIED  TO 
OTHER  SUBSYSTEM  INPUT. 


+  1.00 

+  1.00 

+ 1 . 00 

+  1.00 

♦  1.  16 

+0.  69 

♦  1.00 

+  1.00 

*1.  13 

+0.70 

+  1.00 

♦  1.00 

+  1.  15 

♦0.70 

+  1 . 00 

+ 1 . 00 

♦  1.15 

+0.  70 

+  1 .  17 

♦0.  68 

+  1.  15 

♦0.  70 

+  1.15 

+0.69 

+  1 .  15 

+0.70 

♦  1.  15 

+0.  69 

FOR  THIS  RUN:  REFERENCE  INPUT  *  UNIT  STEP; 

DISTURBANCE  INPUT  *  0.5  STEP  DELAYED 
BY  0.6  SEC  APPLIED  TO  INPUT  QF 
OTHER  SUBSYSTEM  INPUT. 


0 

+0 . 60 

+  1.40 

+2.  30 

♦0.35 

1 

+ 0 . 60 

♦  1.40 

+0.  14 

*  1 . 00 

- 

+0 .  60 

+  1.40 

+0.  34 

+  1.33 

+0 . 60 

+  1.40 

+0.  79 

+ 1 . 2» 

a 

+'■>.  60 

+  1.40 

+  1.14 

+0 . 85 

+«:■ .  6o 

+  1 . 40 

+  1.15 

+0.  70 

* 

♦O .  60 

+ 1 . 40 

+  1.14 

+0.  70 

T 

+0 .  60 

+  1.40 

+  1 .  14 

+0 .  70 

3 

*1.1') 

+0.91 

+  1.14 

+0 . 70 

3 

+  1.16 

+0.67 

+  1.14 

+0.  70 

IV 

*1.15 

+0.  69 

+  1.14 

+0.70 

1  i 

+  1.15 

+0. 69 

+  1.14 

+0.  70 

~QP  rH!S 

pLN:  REFERENCE  INPUT 

=  UNIT  STEP  DELAYED 

0.6  SEC: 

DISTURBANCE  INPUT 

*  0. 5  STEP 

APPLIED 

OTHER  SUBSYSTEM  INPUT. 

'  ’ 

+0 .  6o 

+  1 . 40 

•*■2.  70 

+0.  35 

t 

+  1.17 

+0.81 

*2.  30 

+0.35 

- 

+  1.17 

+0. 66 

+2.  30 

+0.  35 

- 

+  1.16 

♦0. 68 

+3 .  30 

+0.  35 

4 

+  1.16 

+0.  aS 

+  2.  70 

+0.  35 

c 

-J 

+  1.16 

+<1).  68 

♦O.  17 

+ 1  .  oo 

6 

+  1.16 

+0 .  68 

♦0.  42 

+  1.21 

7 

+  1.16 

+o .  cS 

♦0 , 90 

+  1 . 03 

3 

+  1.16 

+'  > .  60 

♦1 .  l£> 

+0. 68 

3 

+  1.16 

+0.  08 

♦l.  15 

+0. 67 

1  V 

+  1.16 

♦< ).  68 

-1 .  13 

+0. 67 

SUBSYSTEM  I 


SUBSYSTEM  2 


FEEDBACK 

SAINS 

FEEDBACK 

GAINS 

SUBSYSTEM  POLES 

UNDAMPED 

NATURAL 

INITIAL  ! 

Kill  : 

FINAL* 

Kll 

K23  : 

K24 

S  • 

FREQUENCY 
< RAD/S) 

TABLE  VII.  DECENTRALIZED  TUNING  EXAMPLE:  TUNING  HISTORIES  FOR 
VARIOUS  TUNINGS 


FOR  THESE 

RUNS: 

SUBSYSTEM  1 

INPUT  -  UNIT 

STEP; 

SUBSYSTEM  2 

INPUT  -  0.5  STEP  DELAYED 

BY  0.6  SEC. 

REFERENCE  RESPONSES:  SECOND-ORDER  STEP 

RESPONSES  (POLES  AT  S< 

—3)  . 

LOCAL  FEEDBACK 

GAINS: 

Kll 

AND  K12 

!  CORRESPOND 

TO 

SUBSYSTEM 

1 

K23 

AND  K24 

'  CORRESPOND 

TO 

SUBSYSTEM 

2 

ITERATION 

I 

Kll 

K12 

K23 

K24 

0 

+0.30 

+  1.30 

+0.25 

+0.  10 

1 

+0 . 87 

♦  1.00 

+0.25 

+0.  10 

-> 

♦o.as 

+0.93 

+0.23 

♦O.  10 

3 

+0.85 

+0.93 

+0.25 

+0.  10 

0 

+2.00 

+0.30 

+0.  30 

0.00 

1 

0.006 

+0.83 

+0.30 

0.  00 

r* 

♦0.03 

+  1.23 

+0.30 

0.00 

3 

♦0.  38 

♦  1.31 

♦0.  30 

0.00 

4 

+0.83 

+0.91 

+0.30 

0.00 

3 

+0.88 

♦0.79 

+0.  30 

0.00 

6 

+0.87 

+0.82 

+0.  30 

0.00 

7 

+0.87 

+0.82 

♦0.30 

0.00 

0 

+2.00 

+0.30 

+0.50 

0.00 

1 

+0.50 

+0.79 

♦O.  50 

0 .  00 

n 

+0.  77 

♦O.  90 

+0.50 

0 . 00 

3 

+  1.03 

+0.76 

+0.  50 

0.  00 

4 

+  1 . 09 

+0.  67 

+0.50 

0 . 00 

5 

+  1.09 

♦0.68 

+0.50 

0 . 00 

+ 1 . 09 

+0.  68 

+0 . 50 

0 . 00 

0 

+0.  50 

+  1.50 

+0.  40 

0 .  00 

1 

+ 1 . 09 

+0.  76 

+0 . 40 

0 . 00 

2 

+ 1 . 04 

♦O .  68 

+0 .  40 

0 . 00 

3 

+ 1 . 04 

+0.69 

+0.  40 

0 . 00 

4 

+  1.04 

+0.69 

+0 .  40 

0.  00 

0 

+0.50 

+  1 . 50 

+0.  40 

+0.  14 

1 

♦  1 . 04 

+0.92 

+0.  40 

+0.  14 

2 

+  1.06 

+0.71 

+0 .  40 

+0 .  14 

3 

+ 1 . 05 

+0 . 73 

+0 .  40 

+0.  14 

4 

+  1.05 

+0.  73 

+0.40 

+0.  1  4 

0 

+0.  50 

+  1 . 30 

+0 .  40 

+0 .  20 

i 

+ 1 . 03 

♦O.  94 

+0 .  40 

+0 . 20 

2 

+  1.07 

+0.  72 

+0 .  40 

+0 . 20 

3 

♦  1 . 06 

+0.74 

+0 .  40 

+0 . 20 

4 

+  1 . 06 

+0.74 

+0 .  40 

+0.  20 
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